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■ ABSTRACT 

. Aims. In this paper we study the feasibility of inferring the magnetic field from polarized multi line spectra using two methods: The 

K ' pseudo line approach and The PCA-ZDI approach. 

Methods. We use multi-line techniques, meaning that all the lines of a stellar spectrum contribute to obtain a polarization signature. 
The use of multiple lines dramatically increases the signal to noise ratio of these polarizations signatures. Using one technique, the 
. pseudo line approach, we construct the pseudo line as the mean profile of all the individual lines. The other technique, the PCA-ZDI 

approach proposed recently by Semel et al. (2006) for the detection of polarized signals, combines Principle Components Analysis 
(PCA) and the Zeeman Doppler Imaging technique (ZDI). This new method has a main advantage: the polarized signature is extracted 
using cross correlations between the stellar spectra and functions containing the polarization properties of each line. These functions 
, are the principal components of a database of synthetic spectra. The synthesis of the spectra of the database are obtained using the 

radiative transfer equations in UTE. The profiles built with the PCA-ZDI technique are denominated Multi-Zeeman-Signatures. 
Results. The construction of the pseudo line as well as the Multi-Zeeman-Signatures is a powerful tool in the study of stellar and solar 
magnetic fields. The information of the physical parameters that governs the line formation is contained in the final polarized profiles. 
^ , In particular, using inversion codes, we have shown that the magnetic field vector can be properly inferred with both approaches 

$_H ' despite the magnetic field regime. 

Conclusions. 
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1. Introduction ization profiles. These mean profiles serves in a secondary 

step as input for codes that produce magnetic stellar surface 

^ The majority of cool stai's, including our Sun, have magnetically j^aps. Observations at different phases of s tar's rotational pe- 

^ confined atmospheres. The study of the magnetic activity is im- riod are required to produce these maps (iBrown et al. 19911: 

. portant for understanding: 1) the dynamo effect, responsible for Ippnati & Brown 1997tlH^n et al. 2000h . 

T-H generating solar and stellar magnetic fields, and 2) the different ^he simplest way to employ the line addition principle is to 

O states of stellai- evolution and the influence of the magnetic field ^^g^^gg spectral lines. However, this brut Une addition tech- 

. during these stages of evolution. nique is not currendy used to retrieve the magnetic field vector, 

7^. The development of the Zeeman Doppler Imaging technique, not even for solar observations, where the inversion of the Stokes 

>■ . described in a series of five papers, is a benchmark in the stel- pai-ameters in individual specti-al lines is a commonly employed 

lar spectropolarimetry domain and, consequendy, in the study of technique. Semel et al. (2009) have recendy shown that the sig- 

magnetic stars by observational methods. It opened a new win- nal obtained from the average of multiple spectral Unes, the so- 

dow in the research of magnetic activity in cool stars through called pseudo line, can be employed to estimate the longitudinal 

the measurement of the circular polarization direcdy from the magnetic field using the center of gravity method (e.g. Rees et 

spectra. The circular and also the linear polarized signals are 1979). 

detectable through a combination of the Doppler and Zeeman :„ this work we present, in Sect. 2, more accurate inversion 
effects ( Semel 1989D. methods applied to the multi line approach by means of the con- 
Few years after the development of the ZDI technique, struction and inversion of the pseudo line, and we extend the 
Semel & Li (1996) proposed the line addition principle to im- study to the linear states of polarization since considering the 
prove the signal-to-noise ratio and make possible the detec- Stokes parameters (Q,U) is necessary to fully determine the ori- 
tion of extremely faint polarization signals in stars. Since then, entation of the magnetic field. 

the addition principle has been used in different and more so- The second part of this work. Sect. 3, is dedicated to an- 

phisticated techniques, and it has been particulai-ly success- other approach of detection of magnetic fields in stai's, where 

ful with the Leas t Squares Deconvolution (LSD) technique all the lines contained in a given spectral range contribute to the 

dDonatietal. 1997h . to construct mean intensity and/or polar- final Stokes profiles. We first present the basis of the so-called 

PCA-ZDI technique. We then describe the procedure to obtain 

Senrf Oj^nTjr re^Me.sfj fo: jramirez@astroscu.unam.mx the final Stokes profiles named Multi-Zeeman-Signatures. We 
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then apply the developed technique to observed spectra of three 
cool stars. Finally, using synthetic spectra, we will show that the 
stellar magnetic fields can be correctly retrieved performing di- 
rect inversions of the Multi-Zeeman-Signature profiles. 
The general conclusions are presented in Sect. 4. 

2. The pseudo line approach 

In this Sect, we concentrate on the line addition technique us- 
ing the simplest way of combine multiple spectral lines that is 
averaging them. We have named the pseudo-line (PL) to the re- 
sulting mean signal because it does not come from any physical 
entity, i.e. it has not associated any intrinsic atomic property like 
Lande factor, potentials of excitation, etcetera. 

The Stokes vector of the pseudo line , Spi - (/, Q, U, V), is 
obtained averaging by separate in each Stokes parameter. 



"I 

z 

1=1 



[li(AR), QijAg), Ui(AR), Vi(AR)] 
ni 



(1) 



where A^ = A - Ac denotes the reduced wavelength and «/ is the 
number of spectral lines. 

In table[T]are listed the individual lines that we have included 
in this work for the construction of the pseudo line. They have 
be chosen to be the same lines as in our previous work Semel et 
al. (2009), hereafter refeiTed as Paper I. 

Considering the solar case of a single point on the stellar 
surface, we used the code diagonal (Lopez Ariste & Semel 1999) 
to compute the Stokes profiles of the lines listed in table [1] and 
consequently, to obtain the profiles of the pseudo line. 

As mentioned in Paper I, we will assume that there is only 
one magnetic field vector on the stellar surface responsible of 
the polarized signals and that the total radiation comes from this 
magnetic. We do not ignore that more sophisticated scenarios 
can be considered, but we are not interested in simulate a real 
stellar spectra but in show that the capabilities of the pseudo line 
to retrieve the magnetic field vector. 

In Fig.[Tl we show two examples of the pseudo profiles com- 
puted using the same given atmospheric model (see Sect. 2.3), 
but using different magnetic strength field. In the left panel, con- 
sidering a field strength of 400G, the profiles of all the individ- 
ual lines are shape similarly, and consequently, the pseudo pro- 
files keep the same shape. In the right panels, considering a field 
strength of 4kG, the similarity in the shape of the individuals 
profiles is not preserved. Note for instance the variety of shape 
in the I and Q profiles. 

The relatively reduced number of individual lines used to 
construct the pseudo line is a first step for a future generalization 
where much more lines can be included. 

Before to continue with the inspection of the pseudo line we 
present the procedure that we follow to express the spectral lines 
in function of the Doppler coordinates. 

2.1. Variable transformation and adequate spectral 
resolution 

In this Sect, we describe a variable transformation from the 
wavelength coordinate into the velocity one: S(A) — > S(X). 

At present, the typical spectrographes used for stellar spec- 
tropolarimetry are cross-dispersed type and cover a range of 
thousands of A. We will employ in our analysis a discrete se- 
lection of lines, but we present the general procedure as if we 
were dealing with the complete spectral range. 



Consider the transformation equation 

X^c logiA/Ao) 



(2) 



where X is expressed in km/s, and A in A. c is the speed of 
light and Aq the minimum wavelength value in a given spectral 
interval. Deriving the precedent Eq., 



dA 

dX = c— , 
A 



(3) 



is possible to obtain a constant step in km/s along the spectra 
when the velocity values are discretized: dX —> AX. 

We have fixed the constant step in AX to I km/s, such that we 
can invert the equation ^ to express A as function of X as: 



A - AoexpiX/c) 



(4) 



where X= [0,1, 2,...]. 

Let Aobs denote the discrete wavelength values obtained in 
the spectra after the data reduction. By imposing a constant step 
in AX one problem arises for certain values because the rela- 
tion A{X) given by Eq. (|4| will not necessarily coincide with the 
values of Aghs- An interpolation is then applied to determine the 
exact value at the wavelength of interest. Let A„ be one of the in- 
exact values such that A„hsj < A„ < Agbsj+i, and letS(X„) denote 
any of the Stokes parameters at wavelength A„. The value S {X„) 
will then be calculated by: 



S(X„)^S(Ai) 



S(Ai^i)-S(Ad 
Ai+i - Ai 



(An - Ai) 



(5) 



With this algorithm, the spectra S (A), observed or synthetic, 
will be transformed to as function of the new variable S (X). 

2.2. The atmospheric model 

In what follows of this Sect, we explore the use of accurate in- 
versions methods to apply to the pseudo line, i.e. inversions tech- 
niques based on the radiative transfer equation (RTE). 

From experience in the solar data analysis, we know that 
some correlation between the atmospheric parameters can gener- 
ate ambiguities in the determination of the magnetic field when 
the noise level is close to the polarized signal levels, (e.g. del 
Toro Iniesta & Ruiz Cobos 1996; Bellot Rubio & Collados 2003; 
Martinez Gonzalez et al. 2006, and references therein). For sim- 
plicity, we initially avoid any possible model parametric degen- 
eration by assuming that the atmospheric model is known and 
only the magnetic field vector is unknown. With this premise. 

Table 1. Central wavelength, transition levels, relative intensities 
and Lande factors of the spectral lines of the multiplet 8 16 of the 
neutral iron. The relative intensities are taken from Allen (2000). 
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Fig. 1. The dashed lines represent the Stokes profiles of each of the individual lines listed in table 1. The solid line represents the 
pseudo line obtained from the rest of the lines. Left side: Stokes profiles considering a moderate field strength of 400 G. Right side: 
The same as in the left panel but considering a strong field strength of 4 kG. In both cases the geometry of the magnetic field is 
azimuthal and inclination angles of 70 and 45 degrees respectively. 



we test the inference of the magnetic field through the inversion 
of the pseudo profiles under an ideal scenario where all the vari- 
ables in the RTE are related to the magnetic field vector : The 
field strength (B), the azimuthal angle {-yazi) the inclination 
field ijinci)- The atmospheric model we chose to compute the 
Stokes profiles is as follows: we fixed the gradient of the source 
function (Ar = 10), the Doppler width to 50 mA, the ratio of the 
line to continuum absorption (77) is twice the relative strength of 
the components of the iron multiplet and the damping parameter 
value is fixed to 0.01. 

The four Stokes parameters for all the individual Fe lines 
are computed considering the same atmospheric model and the 
values of the magnetic field varies randomly in the following 
ranges : B=[0,10] kG, rfl,,=[0,180] degrees and jm = [0,90] 
degrees. 

Consider now a given combination of the magnetic field pa- 
rameters, Y"-(B", y^"^., y,-,"^,). From the synthesis of the Fe lines 
for this combination of parameters, we then obtain the respec- 
tive pseudo line for this Y" model. Note that we are in this way 
associating a magnetic model to the pseudo line. 

2.3. The inversion code 

The exercises consist to determine the magnetic field supposed 
punctual and fixed in a known position in the stellar surface. The 
inversions in each Fe line and in the pseudo line have been done 
separately and the goal is to compare the results. 

We have developed one inversion code per Fe line and one 
for the pseudo line. The code's functionality is : First, we con- 
struct a representative synthetic database of a large number of 
profiles (56600). Then, given a Stokes vector to invert, we find 
the best-fit solution from the database using the same algorithm 
as described in Ramirez Velez et al. (2008). 

We inverted a set of 500 profiles for each of the lines listed 
in Table [U To compute the sets of profiles to invert, we followed 
the same procedure as when constructing the database profiles: 
we have fixed the atmospheric model and we consider random 
magnetic field vectors. 

While inverting the set of Stokes profiles, we applied and 
compared different strategies to recover the magnetic field vec- 



tor Initially, we considered the four Stokes parameters (I,Q,U, V) 
to simultaneously retrieve the solution of the three magnetic pa- 
rameters (B, jazi, Jmci)- We then inverted the same set of pro- 
files using only the I and V Stokes profiles, which is the most 
common case for stellar observations. Finally, we performed the 
inversions using only the two linear Stokes parameters (Q and 
U). 

The goal of performing different inversion strategies is 
twofold. On one hand, we compare the degree of accuracy of 
the inversions using the four Stokes parameters to the most typ- 
ical case in the data analysis that is when we only dispose of the 
intensity and circular polarization spectra. 

On the other hand, this procedure permits to improve the 
precision of the inversions. Let us call the discerning criteria 
as the use of different Stokes parameters to retrieve the differ- 
ent components of the magnetic field. Given that the inversion 
code works finding as solution the best-fit of the profiles in the 
databases, the employment of the discerning criteria has the ad- 
vantage that the databases become largers and thus, the code re- 
duce the error of the inversions. 

We remark however that the results of the analysis of the 
pseudo line do not depend in the proposed criteria described be- 
low, such that some readers may want to skip the next subsec- 
tion. 



2.3.1. Database extension 

To clarify the idea of the extension of the database, let P„ denote 
the Stokes vector for a given set of parameters such that /"„ = 

i'„(y«) and y"=(B«, r;„., ri)- 

Let (B', Jjjj^i) denote the associated solution values found 
with the inversion performed using the I and V Stokes parame- 
ters. Now, let y be the solution value after inversion of the 
and U Stokes parameters of the same P„. Since the values of the 
parameters were retrieved independently, then j + k (in more 
than 90% of the cases). Finally, we use the so-called discerned 
criteria to construct the final solution with a parameters combi- 
nation Y'"''-(B-', y^*j., j/^^i). Since Y""' was not originally in the 
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Fig. 2. First and second columns: Inversions performed using different Stokes parameters (indicated in the Y-axis). Due to the 
inclusion of the magneto-optical effects in the computation of the profiles, the solutions in the second column contain the correct 
information in yi„ci parameters but with wider ranges of errors than in the first column. Third column: All the Stokes parameters are 
used to retrieve the magnetic field values. Last column: Comparison of the error inversions using the discerning criteria (in black) 
and with out it (in gray), more details are given in the text. The better results coiTespond to the discerned inversions. 



database, we are "extending" the database through the use of the 
discerning criteria. 

In Fig.|2] using a set of 500 profiles without any added noise, 
we show the inversion results following the described criteria. 
The graphics correspond to the inversions of the Fe I line at 6141 
A, the first line listed in Table [1] 

While using only the I and V Stokes profiles to perform the 
inversions, first column in Fig. |2] more refined inversions are 
obtained for B and yi„ci that when all the Stokes parameters are 
included in the inversions (third column). This is explained be- 
cause the number of magnetic models whom parameters B and 
jinci are close to the original pair of value parameters {B°", y°"^i) 
get increased when the y,,,, is ignored. 

In the second column we present the results of the inversions 
using only the linear Stokes profiles (Q,U). In this case, also an 
improvement happens for the inversions of the parameter jazi 
compared to the inversions of the third column when the four 
Stokes parameters are considered. 

From the comparison of the inversions, showed in the forth 
column, we then conclude that the discerning criteria improves 
the inversion results: In the case of the field strength, the disper- 
sions in the eiTors diminished considerable and the maximum 
error value decreased from ~140G (symbols in gray) to ~25G 
(symbols in black). In the case of the geometrical parameters. 



the angles y,„c/ and ja-j, the inversions also improve with the 
adopted criteria. 

Note that it is expected that the proposed discerning criteria 
increase the inversion precision only for those codes that em- 
ploys a big set of profiles (databases) where the solution is found. 

2.4. Testing the inversions of the pseudo line 

One advantage of the developed inversion codes is the estimation 
of the typical eiTors. Thus, in order to compare quantitatively the 
inversions retrieved from each line, the results are presented in 
terms of the absolute value of the difference between the original 
and the inferred parameters: 

gori _ gsol ^ (1^51^ |^^^^.|^ \hyi„,l\). (6) 

Additionally, in order to study the effect of the noise in the 
determination of the magnetic field, we add to the profiles three 
different noise levels with respective values of 1x10"^, 1x10"^ 
and 1x10"' of the continuum level (/c). 

In Fig. [3] we show the percentage histograms of the inver- 
sion's eiTors for all the Fe lines (in gray) and for the pseudo line 
(in black). 

In the case without any added noise, we can see that all the 
inverted parameters are well retrieved. The histograms of the er- 
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Fig. 3. Percentage histograms inversions of the Fe Unes (in gray color) and of the pseudo line (in black color). N represents the 
percentage number of inversion cases. In the row direction, ordered by parameter of the magnetic field configuration: the strength, 
the azimuthal and, inclination angles. In the column direction, the comparison of the inversion by noise level. 



rors have a sharp peak at and drop quickly. Note for instance 
that the field strength is inferred with an error smaller than 10 
G with a high probability. The small errors for all the param- 
eters are mainly due to the finiteness of the database, and are 
identified as the precision errors of the database. In any case, 
note that the results with the pseudo fine are the very same the 
ones obtained from the inversions of individual spectral lines. 
We then conclude that the pseudo line encodes the information 
of the magnetic field vector as the individual spectral lines do!. 

In the lowest noise case (10"^/c), all parameters are also cor- 
rectly retrieved, the errors being very close to the precision of 
the database. However, when increasing the noise levels to 10"^ 
and to 10"'/^, the error bars become more important, making 
the error distributions shghtly wider. We remark that for these 
noisy cases the pseudo line is the one that gives the best results, 
corroborating that the addition of multiple Unes increases the 
signal-to-noise (S/N) ratio. 

2.5. Inversions with an uni<nown atmospheric model 

We continue our inspection of the pseudo line relaxing the con- 
straint that the magnetic field is the only free variable. We have 
thus included some line formation parameters in the atmospheric 
model, to verify that in a more general case the goodness of the 
approach remains. 

Considering as free variables the Doppler width (Vd), the 
absorption line ratio (rjo), the source function gradient (Vt) and, 
the magnetic field vector, we have repeated the same exercise as 
before inverting a set of 500 pseudo profiles. 

The construction of the pseudo line follows the same princi- 
ple as before: Given a combination of the atmospheric parame- 



ters (excepting tjq) those values are used to compute all the Fe 
Unes, and subsequently, obtain the pseudo Une. In the case of the 
rjo parameter, a factor (denoted rjj) modifies in the same amount 
any of the t/o values, preserving the relative intensities (i.e. for 
all Fe lines, t/o = t/o • ?//)• 

In Fig. 4, we show the obtained results. After inversion of 
the 500 pseudo profiles, we have found that the retrieved values 
in any of the parameter are correct (with of course and asso- 
ciated error). We then conclude that the atmospheric model and 
the magnetic field can both be deduced from the inversions of the 
pseudo Une. We argue in favour of the results, that the method 
used to invert the pseudo line is the very same the one used for 
the construction. 

With the present results, we then conclude that the multi line 
approach based in the direct addition of many spectral lines, i.e. 
in the construction of the pseudo Une, is a very useful tool in the 
study of solar and stellar magnetic fields. 

3. The PCA-ZDI approach 

Alternatively to the method presented in the last Sect., instead 
of make more efforts in this direction considering as many lines 
as possible in the pseudo line, we prefer to employ a more so- 
phisticated approach recently proposed by Semel et al. (2006), 
hereafter referred as Paper II. In this Sect, we focus in this sec- 
ond approach named PCA-ZDI. 

The principal idea of the PCA-ZDI technique is to incor- 
porate the computation of the radiative transfer equations, for 
instance using an atmospheric model in Local Thermodynamic 
Equilibrium (LTE), in combination with a powerful statistical 
tool, the principal components analysis (PCA), in order to re- 
trieve reliable detections of the stellar magnetic fields. This is 
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Fig. 4. Inversions of a set of 500 pseudo lines performed using the four Stokes parameters. The title in each graphic indicates the 
atmospheric parameters. 



in fact a way to overcome the employment of the weak field 
approximation used originally in the formulation of the ZDI 
technique (Semel 1989) and preserved in the LSD technique 
dPonati et al. 1997h . 

3.1. Detection of stellar magnetic fields with PCA-ZDI 

As will be shown, the PCA-ZDI technique is a powerful 
method for detection of magnetic fields in stars. In a sum- 
marized description of the developed technique, firstly is re- 
quired to dispose of a representative database composed of 
stellar spectra. In fact, despite the domain of application, the 
construction of the database is a key step when is intended 
to perform analysis by pr incipal components ( Rees et al. 2000t 
ISocas Navarro et al. 200 ih. In any case, once defined the 
database, of what ever number of spectra is constructed this one, 
it is decomposed in a new basis using the single value algorithm 
(iGolub & Van Loan 19961) . The new basis is spanned in terms of 
eigenvectors which present many interesting properties. For the 
purposes of this work, we employ the eigenvectors of the new 
basis as detectors to retrieve from the spectra the magnetic signa- 
ture of the Stokes parameters. We have denominated these mag- 
netic signatures as Multi-Zeeman-Signature (MZS). The PCA- 
ZDI technique will thus retrieve a not null MZS profile in the 
polarized Stokes parameters whenever is present a magnetic field 
in the stellar atmosphere and the signal to noise ratio allows it. 

From an observational point of view, the main constrain to 
retrieve the polarized Stokes parameters in individual spectral 
lines is the very faint levels expected in the signals 

3.2. Magnetic detectors 

Based in th e ideas proposed in the first paper of ZDI 
(ISemel 1989h . we have considered that there is only one mag- 
netic element in the stellar surface, at the projected position (fi), 
and we have calculated the local Stokes profiles at this position. 
From the results already obtained with synthetic and observed 
spectra, we advance that the employment of the local profiles 



does not limit the detection in stars with complex magnetic fields 
configurations -dipolar, multi-polar and/or multispot-. A more 
detailed inspection of this statement will be presented in a forth- 
coming paper. 

The stellar spectra in the four Stokes parameters have been 
established with the help of cossaaQ described in Stift (2000) 
and in Wade et al. (2001). cossam is an LTE line synthesis code in 
polarised light that calculates the Stokes profiles over arbitrarily 
large wavelength intervals, both for the sun and for dipolar field 
geometries in magnetic stars. Direct opacity sampling is caiTied 
out over the cr_, the cr+, and the n components separately of the 
(anomalous) Zeeman patterns of the individual lines. 

While covering different positions in the stellar surface and 
at each one considering different field strengths and orientations 
is possible to obtain a representative combination of synthetic 
spectra which will serve to construct the database. 

The position of the stellar element is specified through the 
fi-cos{6Los), where 6los is the angle between the line-of-sight 
direction and the normal to surface element, fi varies from [0.2, 
1] with a step size of 0.2. At each fi position, we have spanned 
the magnetic field strength B=[0,250,...,3000] G and the incli- 
nation angle 7,„c;=[0,15,...,90] °. For economy in CPU-time we 
have fixed the azimuthal angle (jazd^ but note that this fact does 
not limit the detection in the linear Stokes parameters since ap- 
plying rotations of the reference frame to the calculated spectra 
is possible to pass from the Stokes Q to the Stokes U, or to obtain 
any desired configuration in the linear parameters, (e.g. Sect. 6.4 
in del Toro Iniesta 2003). Consequently is justified to consider as 
the general case the approach used in the y,,-, parameter. In this 
way, the combination of parameters produces a database (M) 
with 425 spectra, per Stokes parameter. 

The spectral range covered from 450 to 750 nm, in steps of 
10 mA, producing spectra of 300,000 wavelength points in each 
of the Stokes parameter. In the following, we consider the Stokes 
spectra as vectors S of dimension [300000]. 
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Finally, three different atmospheric models have been con- 
sidered T=3500, 4750 and 5750 °K. For each model we con- 
struct the respective database. In other words, T is fixed to the 
value To in each database. In all cases, we have considered the 
same atomic solar abundances as in Grevesse & Sauval (1998), 
leaving out the molecular solar lines. 

Let represent a given combination of the model parame- 
ters 

= ^,n(B,yinci,yazhl^,To),„, m = [0, 1, ..,424] (7) 

and let S,„(^,„) denotes the mth combination in any of the Stokes 
parameters calculated with cossam. Grouping all the 425 con- 
sidered spectra models, we obtain an array M of dimensions 
[300000, 425] per Stokes pai-ameter 

Applying the Single V alue Decomposition (SVD) procedure 
(iGolub & Van Loan 19961) . the new basis with respective set of 
eigenvectors [P] = [Pq, P\,..., Pa2a\ that span the database M is 
found. 

Consequently, any spectra Sm can be expressed as a linear 
combination of the eigenvectors. Considering the change of vari- 
able previously described in Sect. 2. 1, the correspondent expres- 
sion for the spectra Sm(X) in terms of the eigenvectors is : 

424 

5„,,,(^,„, x)^Yj <i -P".' W ' = Q' ^' (8) 

n=() 

where X expressed in km/s has replaced the wavelength value. 
Note that any Stokes spectra Sm can be already identified by an 
unique associated set of coefficients, {a"']-[a'Q, a'", a"24^- 

In figure |5] we show a small portion of the first three eigen- 
vectors of the database constructed for T=4750 ° K. 

3.3. Space dimension reduction 

One of the reasons why we employ PCA is because of its capa- 
bility of data compression and space dimensions reduction. We 
mean that in Eq. ([8]l the equality is reached only when n^ax - 
424, but in practice «„„, ^ could be dramatically reduced such that 
only the first components are considered. The consequent ques- 
tion is how to determine the number of components really sig- 
nificants i.e. how to determine n„ax- To answer this question we 
have reconstructed all the profiles with 1 component, then with 
2, and so on, and we compare the reconstructed spectra with the 
original ones. The mean percentage contribution per component 
is plotted in Fig.|6] 

From Fig. |6] we can appreciate that rate of contribution of 
the components to the original spectra decreases faster than an 
exponential, permitting to truncate the expansion in Eq. (|8]l to 
the first components. Given that the contribution for n = 9 is 
already inferior to 0.001% in the intensity Stokes parameter and 
inferior to 0.05% in the polarized ones, we have decided to keep 
only the first ten terms in Eq. ([8]l. The same number of compo- 
nents will be used when obtaining the Multi-Zeeman-Signatures 
and when inverting them. Moreover, for any star with the same 
temperature that the one considered in the database, the same ten 
eigenvectors will be useful to retrieved the MZS. 

At this point in the development of the technique we have 
followed a typical approach with PCA. In the following we em- 
ploy the eigenvectors as magnetic detectors such that combining 
them with the ZDI principles, the magnetic signature from the 
spectra will be extracted. 
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3.4. Tine Multi-Zeeman-Signatures (MZS) 

Another of the advantages of the employment of PCA apart from 
the data compression, is that it permits to retrieve a set of coef- 
ficients {a'" ] associated univocally to each one of the spectra 
5m. We will use then these coefficients to establish a relation to 
obtain the Multi-Zeeman-Signatures. 

Given that the eigenvectors have the property of being or- 
thonormals, from Eq. ([Sjl it is then trivial to obtain the set of 
coefficients for any of the synthetic Stokes spectra as: 

a';ii^S„M,„,X)- PnjiX); (9) 

where n= [0,1,... ,424] and i=(I,Q,U,V). The relation obtained in 
Eq. (|9]l represents the fundamental principle that will be em- 
ployed to extract the magnetic signature from the observed spec- 
tra (S"* ')- In order to do it, the last step of the procedure is to con- 
sider the rotation velocity of the star and the associated Doppler 
effect, i.e. combine the MZS profiles with the ZDI principles. 

It is important to mention that in order to equalize the spec- 
tral resolution of the observations with the one of the eigenvec- 
tors, we degraded the resolution of the eigenvectors to the spec- 
tral resolution of the observations. This last is done through a 
linear intrapolation similar the one described in Sect. 2.1. 

We proceed now to apply the same operation expressed in 
Eq. (|9]l to the observed spectra 5"**, but when the eigenvectors 
are Doppler shifted by an amount Y-AX : 

a''„'i\Y) ^ Sf\X) ■ P„j(X-Y). (10) 

The vectorial product in the right side of in Eq. ([TOl i that 
permits to retrieve the or-coefficients at each considered F-value 
is equivalent to a cross correlation function between the spectra 
and the magnetic detectors (the eigenvectors). 

On the other hand, it is convenient to mention that the ex- 
pression found in Eq. ( [TOl i can be replaced by another one where 
the eigenvectors (P) are substituted by the absolute value of the 
same eigenvectors, denoted by { | /" | } : 

<7(F) = Sf"(X) ■ \P„j\(X-Y). (11) 

Considering the whole range of Doppler velocities, 
Y-[Yo,Yi, ...,Ymax], the general expression of the Multi- 
Zeeman-Signature profiles (Pmzs) is given by: 

iPMZs)nAY) = «r(>'0), <.^(Fl), ... , ^'d-,™.-)) . (12) 

Note that for each Stokes parameter, (I,Q,U,V), 425 Multi- 
Zeeman-Signatures are retrieved. Since the eigenvectors are or- 
thonormals, then each one of the 425 Multi-Zeeman-Signatures 
are in principle independents, making possible to perform inde- 
pendent analysis with each one of them by separate. 

3.5. Magnetic detections 

Once presented the technique, we will now apply it to the real 
observed spectra. We will employ the eigenvectors obtained for 
the case of single points in the stellar surface, however this fact 
does not limit the detection of the global magnetic field. 

We will first compare the MZS obtained from Eq. 10 to those 
retrieved following Eq. 1 1 . Let v,„ax and v„„„ denote the maxi- 
mum and minimum projected rotation velocities of the star vsini, 
such that the star velocities span as V"' - [v„„„, v„,ax]- If the 
considered Doppler velocity Y does not coincide with any of the 
V"' values, the intensity MZS profile recovers the continuum 
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Fig, 5. Plots of the first three eigenvectors, in the four Stokes parameters, in a small spectral range of /l=[4500, 4503] A. 
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Fig. 6. Percentage contribution of the first 10 components from Eq. ([8]) in logarithmic scale. Note that the contribution for the 10th 
component in the reconstruction of the original spectra is already at least inferior to 0.05 percent. 
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Fig. 7. Comparison of the MZS for the star Ilpeg considering the absolute value of the eigenvectors (upper panels) and without it 
(lower panels). 
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Fig. 8. Examples of normalized intensity and circular polarization profiles of two solar type stars obtained with the presented PCA- 
ZDl technique. The left plots shows the results for the HD 155555 binary system observed at the the AAT telescope (Australia) in 
April 2007. The right plots correspond to the results for the HD 190771 star observed at the Bernard Lyot telescope in August 2007. 



value, and the polarized MZS profile an aleatory sequence of 
values around zero. 

Contrary, if Y coincides with one of the V"' values, the mag- 
netic signature in the polarized parameters and the total intensity 
profile are retrieved. In Fig. |7]we compare the first nine MZS 
considering the absolute value of the eigenvectors (upper pan- 
els) to the first nine MZS when the absolute value is not included 
(lower panels). The data correspond to the Ilpeg star observed in 
August 2007 with the Bernard Lyot telescope, at the Pic du Midi 
Observatory, using the NARVAL polarimeter 

From Fig. |7] we appreciate that when the absolute value is 
included the nine polarized MSZ are clearly detected with a high 
S/N ratio. If the absolute value is not considered the level of the 
first polarized MZS is clearly superior to the noise level but the 
level of the MZS began to decrease for the rest of eigenvectors. 

A discussion about how to take advantage of this multiple 
Multi-Zeeman-Signature detections will be presented in a forth- 
coming paper In the following, we prefer to add the nine indi- 
vidual Multi -Zeeman-Signatures in order to obtain a total profile 
per Stokes parameter. 

In the last subsection we have showed that the expansion 
in Eq. (O can be reduced to the first ten components. We then 
add the contribution from these components to find the final ex- 
pression that will be employed to retrieve the Multi-Zeeman- 
Signature: 

9 

(PmzsUY) = 2 {<j\Yo), <^.-'(yi), ... , <.^(F„«.)) . (13) 
«=o 

Considering the absolute value of the eigenvectors and fol- 
lowing the procedure of the last equation, in Fig. [8] we show 
two examples of the MZS profiles obtained in intensity and 
in circular polarization of two cool stars. The binary system 
HD 155555, left panels, was observed in April 2007 with the 
Anglo Australian Telescope, at the AAO observatory, using the 
SEMPOL polarimeter The solar type star HD 1 9077 1 , right pan- 
els, was observed in August 2007 with the Bernard Lyot tele- 
scope, at the Pic du Midi Observatory, using the NARVAL po- 
larimeter 

The detection of magnetic fields in this three stars have 
been also found through the LSD method. In fact, the MZS 



profiles from Fig. [8] seems similars in shape to those retrieved 
with the LSD method (Dunstone et al. 2008, Petit et al. 2008). 
Nevertheless a proper comparison has not been done to dis- 
tinguish possible differences in the profiles and how important 
could be those differences when retrieving the magnetic field. 

In any case, for the moment, the PCA-ZDI technique has 
been applied with success to retrieve mostly the circular states of 
polarization in cool stars (Ramirez Velez et al. 2006). In fact this 
is in part due to the absence of representative samples of linear 
data to analyse, so we leave for a future work studies where the 
linear states of polarization can be incorporated in the detection 
and measurement of stellar magnetic fields. 

3.6. Inversions of the Multi-Zeeman-Signatures 

In this section we intend to show that it is possible to infer the 
magnetic field through the inversion of the MZS profiles. We will 
place a simplistic scenario where again we consider local Stokes 
profiles from different surface elements. We further assume that 
there is a punctual magnetic field in the surface elements. From 
now on, the MZS to which we refer will be obtained from the 
local Stokes profiles through Eq. (13). Given a MZS, the inver- 
sion exercise consist in find the position of the surface element 
and the magnetic field vector. 

The local Stokes profiles are calculated for a small spec- 
tral range [5400,5500] A. The parameters of the Stokes profiles, 
and thus those of the MZS, varies randomly in the following 
ranges : the magnetic field strength [0,10] kG, the inclination 
field [0,90]°, the magnetic azimuthal angle [0,180]° and the pro- 
jected position of the surface element (yu) [0,1]. 

Since the inversions are performed in the space of the MZS 
we have produced 6000 MZS that serves as database where we 
tested the inversion of a set of 600 MZS. 

As mentioned previously there are two alternatives to pro- 
duce MZS, namely, whether or not is considered the the abso- 
lute value of the eigenvectors (Eqs. 10 and 1 1). We will test both 
alternatives and we will show that whether or not is included the 
absolute value, the inversions are correct. 

We first present the case where the absolute value of the 
eigenvectors is employed to retrieve the MZS. In upper panel of 
Fig.|9]we show an example of a synthetic MZS (in black color) 
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and the respective MZS solution (in gray color). The inversion 
is correct not only because the similarity between the MZSs but 
because the values of the parameters solution are close the orig- 
inals ones (see the title in each panel). Moreover, given that any 
MZS is associated to a stellar spectrum, in the lower panels of 
Fig.|9]we corroborate that, as expected, the Stokes profiles asso- 
ciated to the MZS solution fit enough well the originals ones. 
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Fig. 9. In black the original profiles, in gray the retrieved so- 
lution. Upper panels, examples of the MZS when the absolute 
value of the eigenvectors is considered. Lower panels, a small 
spectral range of 10 A of the Stokes profiles associated to the 
MZSs. 



In Fig.[TO]we show the MZS associated to the same combina- 
tion of parameters as before but for the case where the absolute 
value is not considered. The solution parameters in this case are 
also close to the original ones (upper panels) and the associated 
solution Stokes profiles gives a good fit to the original Stokes 
profiles (lower panels). 

Finally, we present the results of the inversions of the 600 
MZS. In Fig. [TT] the upper panels correspond the case when 
the absolute value of the eigenvectors is included to produce the 
MZS and the lower panels correspond to the case when is not 
included. 

These results show that in both cases, the magnetic field vec- 
tor is correctly retrieved despite the position in the stellar sur- 
face, for all strength fields and orientations. 



The final conclusion in this section is that the Multi-Zeeman- 
Signature represents not only an effective tool for the detection 
of stellar magnetic fields, but also that inversions in the space of 
the Multi-Zeeman-Signatures can be directly performed. 



4. Conclusions 

In this article, we have focused in the multi line analysis of spec- 
tropolarimetric data, with emphasis in the detection and infer- 
ence of stellar magnetic fields. The use of multiple lines is re- 
quired to increase the signal to noise ratio in those cases where 
the polarized signals in individual spectral lines are significantly 
inferiors to the noise level, remaining hidden for analysis pur- 
poses. This is typically the case, for instance, in observations of 
cool Stars. 

To deal with this observational constraint, we have presented 
the development of two new approaches that serve both of them 
to extract a polarized signal from a mixture of several individual 
spectral lines. While employing any of the proposed multi-line 
techniques, namely the pseudo line or the PCA-ZDI approach, 
we remark three main advantages: (1) the contributions of all 
the spectral lines to the final polarization signal are included - 
despite the similarity or not in the individual shape profiles-, (2) 
it is a valid method for the detection of the circular and the linear 
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Fig. 10. Same as in Fig. |9] but when the absolute values of the 
eigenvectors is not considered to produce the MZS. 
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Fig. 11. Inversions results of the entire set of MZS profiles. The 
upper (lower) panels are the results for the case when the abso- 
lute value of the eigenvectors is (not) included to produce the 
MZS. In both cases the parameters are properly retrieved. 



states of polarization and (3) it is not limited to weak or to strong 
magnetic field regimes. 

In the case of the pseudo line, we have presented an inversion 
code that builds a database of pseudo lines each one attached to 
a particular atmospheric model. 

Initially, considering an ideal scenario where the magnetic 
field is the only free variable of the atmospheric model, we have 
found that the results obtained from the inversions of the pseudo 
line are as goods as those obtained from individual spectral lines. 
We have also shown that it is expected that in the case of real 
observed spectra (with a given noise level), the best results in 
any of the magnetic components are retrieved with the inversions 
of the pseudo line. Finally, we have showed that the atmospheric 
model, in addition to the magnetic field, can also be inferred 
through the inversions of the pseudo line. 

In the case of the Multi-Zeeman-Signatures, we discussed in 
detail the development of the technique, and we applied it to real 
observed data in order to illustrate the type of profile obtained 
with the PCA-ZDI technique. 

The circular polarized MZS profiles, like those in Fig. [8] 
where the polarized signal level is clearly superior to the noise 



level, represent in a unambiguously way a detection of the mag- 
netic field present in the stellar object. 

By construction, the Multi-Zeeman-Signatures contain the 
information of the magnetic field and atmospheric model. To 
verify this, we have performed inversions in a synthetic set of 
MZS profiles. The results obtained after the inversions show that 
the all the considered parameters in the atmospheric model are 
correctly retrieved. 

It is pertinent to mention that some refinements have to be 
done to the technique in order to achieve magnetic fields mea- 
surements through the inversion of the MZS profiles. In particu- 
lar, the Doppler broadening effect due to the rotation of the stars 
and the case of inversions of continuum fields distributions over 
the stellar surface will be presented in forthcoming papers. For 
the moment, we have settled the basis of the approach founding 
that PCA-ZDI is a robust technique for the analysis of stellar 
magnetic magnetic fields. 

The final conclusion of this work is that when the inver- 
sion of the magnetic polarized signal are done through the same 
method used to construct them, the physical parameters that de- 
termine the line formation process can be properly recovered, 
included in particular, but only, the stellar magnetic field. 
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